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Abstract 

Finite time blow-up is shown to occur for radially symmetric solutions to a critical quasilinear Smoluchowski- 
Poisson system provided that the mass of the initial condition exceeds an explicit threshold. In the su- 
percritical case, blow-up is shown to take place for any positive mass. The proof relies on a novel identity 
of virial type. 



^ ' 1 Introduction 

■ We study the occurrence of blow-up in finite time for radially symmetric solutions to a generalized Smoluchowski- 
^ Poisson system 

d 



OO 

o 



d t u = div (a(u) S/u - u Vv) in (0, oo) x 5(0, 1), (1) 
= Av + u-M in (0,oo) x 5(0, 1), (2) 
d u u = d u v = on (0,oo) x 95(0,1), (3) 

u(0) = u Q > in 5(0,1), / v(t, x) dx = for any t G (0, oo), (4) 

Jb(o,i) 

where 5(0, 1) denotes the unit ball of M. n , n > 2, and M the mean value of uq. The diffusion coefficient a 
belongs to C 2 ([0, oo)) and is assumed to be positive for simplicity (see Remark [5] below) . The system (P)-(jl|) 
arises in astrophysics as a model of self-gravitating Langevin particles jl] and in biology [5] where it is also 
known as the parabolic-elliptic Keller-Segel chemotaxis model. 
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A fundamental property of solutions to (|TJ)— d3J) is that 



u(t)>0, / u(t, x) dx = M |.B(0, 1)| , and 

JB(0,1) 



5(0,1) 



v(t, x) dx = 



(5) 



for £ > 0, which readily follows from (pQ), ([3|), the comparison principle, the non- negativity of uq, and the 
definition of v. It is by now well-known that, if a = 1, there are radially symmetric initial data uq for which 
the corresponding solution to (HI)-© blows up in finite time and this singular phenomenon may happen for 
any M > if n > 3 but only for M > 8ir if n = 2 [HE]. More recently, it was shown in [3] that there is 
a critical exponent for the nonlinear diffusion coefficient a which separates two different behaviours for the 
solutions to dH)-(jl]): on the one hand, if a(z) > c {l + z) a and a > (n — 2)/n, there is a unique global classical 
solution to (HI)-© for any non-negative initial condition no € L°°(B(0, 1)) (and this is actually true for a 
general smooth bounded domain of R n , n > 1). On the other hand, if a(z) < c (1 + z) a and a < (n — 2)/n, 
radially symmetric solutions to (pQ)-(|4|) blowing up in finite time are constructed in [3]. 

Except for n = 2, the critical case a = (n — 2)/n is not covered in [3j and the purpose of this note is to 
fill this gap: indeed, the main outcome of our analysis is that, if a(z) < c (1 + z)( n ~ 2 ^ n , there are solutions 
to ([I])-© blowing up in finite time when M exceeds an explicit threshold. As a by-product of our analysis, 
we also establish an alternative and simpler proof of the blow-up result in [3] for a G [0, (n — 2)/n). Indeed, 
the construction of solutions to ([I])-® blowing up in finite time performed in [3] relies on the possibility of 
reducing ©-([l]) to a single parabolic equation. The approach used in this paper is completely different and 
relies on the derivation of a differential inequality of virial type which cannot hold true for all times. When 
a = 1, this technique is used in [6j where it is shown that the moment of order n of u satisfies a differential 
inequality which contradicts the non- negativity of u after a finite time. Seemingly, the moment of order n 
of u does not give valuable information when the diffusion is nonlinear and we introduce nonlinear functions 
of u to be able to handle this case. We finally point out that the above results are only valid for n > 2: the 
situation is qualitatively different in the one dimensional case n = 1 and will be considered in a separate 
paper [2]. 

2 Finite time blow-up 

We first introduce some notation: for (a,p) € [0, oo) x (l,oo), we define 





We also define m p (f) for p > 1 and / G L 



(X 



(0, 1) by 




(7) 
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Our main result then reads as follows. 



Theorem 1 Assume that there are a G [0, (n — 2)/n], and positive real numbers c\ > and c 2 > such 
that 

< a(z) < ci z a + c 2 /or z > . (8) 

Let M > and consider a non-negative radially symmetric function uo G L°°(B(0, 1)) suc/i t/iai ||iio||i = M. 
Assume further that EM,p(m p (uo)) < for some p > 1, where 

/ M \ ( 2 P+™ a (P+l))/n / M \ (2p)/n 

1 / A f \ p ^ 

+ M z - ; — , z > . (9) 

p{p + 1) V " J 

Then the system ([I])-® has a unique maximal classical solution (u,v) with finite maximal existence time 
T m ax G (0, oo) and ||u(i)||oo — ► oo as t -> T ma2 .. 

Remark 2 There are initial data u$ for which EM,p{rhp{uo)) < /or some p > 1. Indeed, observe that 
E M , P {°) <0forallM>0ifae [0,(n-2)/n) and /or ikf > n (c\p(p + l)K p ((n — 2) /n)) n//2 ifa = (n-2)/n. 
It is then sufficient to take uq sufficiently concentrated near x = so that rh„(uo) is close to zero (for instance, 
uo = M5~ n 1_b(o,5) f or <5 > sufficiently small). 

Proof of Theorem Q3 

By Theorem 1.3], there exist a maximal existence time T max G (0,oo] and a unique radially symmetric 
classical solution (u,v) G C([0,T max ); L 2 (B(0, 1); M 2 )) n C 1 ' 2 ((0, T ma;c ) x B(0, 1);M 2 ) to ©-(gD satisfying © 
for £ G [0, T m ax)- Moreover, if T max < oo then ^(t)!!^ — ► oo as t — ► T max . 
We introduce 

If 1 I 1 I ' M \ p 

U ^ r ) := ip/n iM / U ^' X ) dx &nd := - / U & r ) r dr 

n|5(0, 1)| / B (o,0 P Jo \n 

for (t, r) G [0 iT max ) x [0, 1] and derive the following identity of virial type. 

Lemma 3 Let A be defined by A' = a and A(0) = 0. Then 

dm v %e 1 /M\ p+1 m . . 

v - = Mm p - — + 7£ p (u) (10) 



dt p(p +1) \ n J 

with 

p-2 



Kp(u) := jf r 2n " 3 - uj ^2(n - 1) ^ - - (p - 1) r n u^) dr . 
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Proof. Integrating (pQ) gives that U solves 



Consequently, 



d t U = r"- 1 d r A(u) + u U - — r n u with U(t, 0) = U(t, 1) = . 

n n 





'M 




v n 


1 


(M 


p 


\ n 






+ /' 


2(n 


JO 





p-1 



M 



(r n - 1) d r Z7 + ( — - U ) d r U - r 2 ^" 1 ) d r A(u) 



\ p+i 

+ M m, + 



1 /M\ 



p + 1 \ n J 



p+i 



dr 



t4(m) dr . 



and hence ([10 



□ 



The next step is to estimate TZ p (u) in terms of m p . To this end, we notice that the assumption ([8]) on a 
warrants that 



< A(z) < 



n z 1+a + c 2 z, z>0. 



l + a 



(11) 



In view of (jll|) . we have 



2(/> - li ( —-U )-(/>- I) r" » 



< max|2(n-l) - C/^ - (p - 1) r n u, o| (^^^ + c 2 



< max <j 2(;> - 1 ) ( Z7 ]—(/>— 1. ) /■" </.() 



ci 



l + a 



2(n-l) / M 



< 



ci(p-l) / 2(n-l) 
l + a V P ~ 1 



l+a 



n 



l+a 



(p — l)r n \ n 
M 



U 



+ c 2 ] u 



u + 2{n-l)c 2 ( — -U) u, 



and thus 



r 1 f m \ p_1 

+ 2(n-l)<52 J r n - 2 I — -U) 0,V<h. 



(12) 



Since a € [0, (n — 2)/n], the function r i— > r^ n 2 an )/ n is concave and we infer from the Jensen inequality 
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(with measure ((M/n) - U) p+a - 1 d r U dr) that 

r n-2-an | _ _ ( r j j- (]/ . 







p + a \n J I Ip + a \n 

< (np )(n-2-g,)/n , M y2p + an(p + l))/n 

p + a \ n J P 

Arguing in a similar way to estimate the second integral in the right-hand side of (fT2j) , we deduce from ([12 
that 

1 (M\ p+1 
n p {u) < E M>p (m p ) -M m p + p{j) + ^ ( — ) 

Inserting this estimate in fjlOH we arrive at 



dm p (t) 
eft 



< EM, P (m p (t)) for i G [0, 

) • (13) 



Assume now for contradiction that Trnax 

= oo. Since z i— > Em, p (z) is an increasing function and 
m p (0) = fh p (uo), we realize that, as soon as EM, P (m p (0)) < 0, we have m p (to) = for some to G (0, oo). 
Thus, C/(t0) r ) = M/n for all r G [0, 1] which contradicts the fact that £/(ioi 0) = 0. Consequently, T max < oo 
and the proof is complete. □ 

Remark 4 If a = 1 (i.e. a = Lemma [5] is a/so wa/ic? /or p = 1 and mi coincides with the moment of 
order n used in |6j/. Not surprisingly, if n = 2 and a = 1, we have Em, p (0) < if M > 4(p + 1) which 
converges to 8 as p approaches 1 and we recover the well-known threshold condition ||«o||i > 87r for finite 
time blow-up to occur in the parabolic- elliptic Keller-Segel system fdjj. 

Remark 5 The requirement a > is only used to have classical solutions to (H])-(jl]) but does not play any 
role in the blow-up condition and the identity of virial type (Lemma\3^. Thus, Theorem^ remains valid if 
the diffusion is degenerate (for instance, a{r) = mr m ~ x with m G [l,2(n — l)/n]) provided an appropriate 
notion of weak solutions is available. 
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